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Abstract: 

We use the quantum action to study quantum chaos at finite temperature. We present a 
numerical study of a classically chaotic 2-D Hamiltonian system - harmonic oscillators with 
anharmonic coupling. We construct the quantum action non-perturbatively and find tem- 
perature dependent quantum corrections in the action parameters. We compare Poincare 
sections of the quantum action at finite temperature with those of the classical action. 
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1 Introduction 



Deterministic chaos has been observed in a huge number of phenomena in macroscopic i.e. 
classical physics. But chaotic phenomena were also found in microscopic physics ruled by 
quantum mechanics - an area presently under very active investigation. For reviews see 
Refs.JT], 0, H|, U ^|. For example, the hydrogen atom in a strong magnetic field shows strong 
irregularities in its spectrum ||. Irregular patterns have been found in the wave functions 
of the quantum mechanical stadium billard [J/]]. Billard like boundary conditions have been 
realized experimentally in mesoscopic quantum systems, like quantum dots and quantum 
corrals, formed by atoms in semi-conductors ||. 

Classical chaos theory can not simply be taken over to quantum physics (due to Heisen- 
berg's uncertainty relation). Hence workers have tried to characterize quantum chaos in 
ways alternative to classical chaos. A fruitful approach is Gutzwiller's periodic orbit theory 
[l] which predicts level densities in the semi-classical regime (Rydberg states in atoms). 
Wintgen || applied it to the diamagnetic hydrogen system and was able to extract periodic 
orbit information from experimental level densities. Another successful route is the analysis 
of level densities in the context of random matrix theory. A conjecture by Bohigas et al. 
states that the signature of a classical chaotic system is described by random matrix 
theory and its spectral density follows a Wigner distribution. 

Strictly speaking, there is no universally accepted definition of quantum chaos. Many 
roads have been travelled searching for "irregular" signals in quantum systems. In the quest 
for quantum chaos, it has been popular to examine systems which exhibit classical chaos. 
However, Jona-Lasinio et al. flHj found chaotic behavior in a quantum system without 
classical counterpart, i.e., a quantum many-body system undergoing multiple resonant 
tunneling. Jona-Lasinio and Presilla [|TTJ further studied quantum many-body systems 



of interacting bosons in the thermodynamic limit. They found chaotic signals when the 
interaction strength drives the system away from integrability. 

Signals of quantum chaos have been investigated also in fermionic many-body systems, 
in particular by looking at the Poisson versus Wigner level spacing distributions |T2], [TB| ]. 
Georgeot and Shepelyansky 0, [15], [31| studied a spinglass model of spin 1/2 particles 
and looked at the transition from the integrable to the chaotic regime, when tuning a 



random transverse magnetic field. Casetti et al. [|T^] considered the large N limit of N- 
component 4 oscillators in the presence of an external field. They compare classical chaos 
with quantum chaos, computed in mean field theory. They observe a strong suppression of 
quantum chaos due to quantum corrections moving the system away from a hyperbolic fixed 



point responsable for classical chaos. Habib et al. [18] have studied in chaotic systems the 
relation between classical and quantum expectation values. They found divergence between 
the two after a time scale t% ~ ln(C/h). Decoherence reduced this discrepancy. 

However, the study of quantum chaos lacks from the following shortcomings: (i) In 
quantum chaos, we have no "local" information of the degree of chaoticity, being available 
in classical chaos via Lyapunov exponents and Poincare sections from phase space, (ii) Also 
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little is known about the role of temperature in quantum chaos. For example, the analysis 
of level densities is insensitive to temperature. To overcome those problems, we adopt the 
strategy of building a bridge between quantum physics and classical physics. This means 
a relation involving the quantum transition amplitude and the classical action. Starting 
from the path integral, the following "bridges" have been considered: 

(i) Sum over classical paths. The quantum mechanical transition amplitude can be 
expressed via the path integral. In certain cases this path integral can be expressed as a 
sum over classical paths only. This holds, e.g. for the harmonic oscillator. Unfortunately, 
such relation holds only in a few exceptional cases [19 |. 

(ii) Gutzwiller's trace formula. Gutzwiller [|J has established a relation between the 
density of states of the quantum system and a sum over classical periodic orbits (periodic 
orbit quantisation). The trace formula has been applied in the semi-classical regime (e.g. 
highly excited states of atom). Wintgen || applied it to the diamagnetic hydrogen system 
and extracted periodic orbit information from experimental level densities. 



(iii) Effective action. Another possibility is to use the convential effective action |20| , [21 



An effective action has been also considered at finite temperature El . The effective action 



has been introduced in such a way that it gives an expectation value < 4> >= 4>dass which 
corresponds to the classical trajectory and which minimizes the potential energy (effective 
potential). Thus one can obtain the ground state energy of the quantum system from its 
effective potential. Because the effective action has a mathematical structure similar to 
the classical action, it looks like the ideal way to bridge the gap from quantum to classical 
physics and eventually solve the quantum chaos and quantum instanton problem. Cametti 
et al. |23| have computed in Q.M. the effective action perturbatively via loop (%) expansion 
and studied quantum chaos in the 2-D anharmonic oscillator. Perturbation theory forces 
to terminate such series at a finite (low) order. The difficulty lies in the estimation of the 
remainder. This poses a problem in the context of chaos, because chaotic systems are known 
to exhibit sensitivity not only to initial conditions but also to parameters of the system. 
There are higher loop corrections to the effective potential V e ^ as well as to the mass 
renormalisation Z. The kinetic term of the effective action is given by an asymptotic series 
of higher time derivatives (infinite series of increasing order). The problem in interpreting 
T as effective action is that those higher time derivatives require more intial/boundary 
conditions than the classical action. 

(iv) Bohm's interpretation of quantum mechanics. An interesting approach has been 
taken by Iacomelli and Pettini [23|. They have used Bohm's |2S| interpretation of quantum 
mechanics, which expresses the Schrodinger equation in polar form. The radial part satisfies 
a continuity equation of a velocity field and the phase satisfies a Hamilton- Jacobi equation. 
The "trajectories" of the quantum system are equivalent to Lagrangian trajectories of 
the fluid velocity field. This allows to apply the concept of classical chaos theory and 
in particular compute Lyapunov exponents. In Ref.plfl this framework has been used to 
study the chaotic behavior of the hydrogen atom in an external magnetic field and different 
regimes of quantum chaoticity have been observed. 
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2 Quantum Action 



In the following we will explore the quantum action |26 as bridge between classical and 
quantum physics. This action has a mathematical structure like the classical action, and 
takes into account quantum effects (quantum fluctuations) via renormalized action pa- 
rameters. The quantum action has the following attractive properties [27]]: (i) It can be 
computed non-perturbatively, via numerical simulations, (ii) It is defined for arbitrary 
finite temperature, (iii) In the zero-temperature limit analytic expressions exist for the 
quantum action, provided that the classical potential is sufficiently symmetric (parity sym- 
metric in D=l and rotationally symmetric in D=2,3). (iv) In a numerical study of the 1-D 
anharmonic oscillator, we searched for the presence of higher order time derivative terms 
in the kinetic term of the quantum action. No signal for such term was found | 27| . The 



quantum action has proven useful for rigorous definition and quantitative computation of 

Let us recall the definition of the quantum action, as proposed 



quantum mstantons 
in Ref.|2§ 



28 



Conjecture: For a given classical action S = J dt^x 2 — V(x) there is a quantum action 
S = J dt^x 2 — V(x), which allows to express the transition amplitude by 



G{xf h tfi\x in ,t in ) = ^exp[- S[x d ] 



Xfi,tfi 



j l n j ( i n 



(1) 



Here x c i denotes the classical path, such that the action S(x c i) is minimal (we exclude the 
occurrence of conjugate points or caustics). Z denotes the normalisation factor correspond- 
ing to S. Eq.(|l|) is valid with the same action S for all sets of boundary positions Xfi, x in 
for a given time interval T = tfi — t in . The parameters of the quantum action depend on 
the time T. Any dependence on xji,Xi n enters via the trajectory x c \. Z depends on the 
action parameters and T, but not on Xfi, Xi n . 

The quantum action at finite time T can be interpreted as action at finite temperature. 
The partition function in quantum mechanics requires to go over from real time to imaginary 

we have shown that the 



time and to impose periodic boundary conditions In Ref.[ 
expectation value of a quantum mechanical observable O at thermodymical equilibrium 
can be expressed in terms of the quantum action along its classical trajectory from Xj n , 
(3 in = to x, (3, given by 



S, 



(3[Xel. 



\x,/3 



+ V(x 



cl 



5S 



5x 



8 l X d ] 



l i 



Xci{/3) 



X 



(2) 



where T denotes the absolute value of imaginary time. The parameters (3, temperature r 
and T are related by (3 = = T/%. In numerical studies with the quartic potential in 1-D 
we found the parameters of the quantum action to be temperature dependent [06[ [27[ . 
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3 Analytic form of quantum action at zero tempera- 
ture 



Why are we interested in the zero-temperature limit? First, from the experimental point 
of view, the ground state properties of a quantum system are often much better known 
than excited states. Second, zero temperature corresponds to large time and from the 
point of view of chaos theory, it is interesting to study dynamical behavior at large times. 
Examples: (a) Behavior of particles moving in storage rings, (b) In classical Hamiltonian 
chaos, there are quasi-periodic trajectories, remaining for a long time in the vicinity of 
regular islands. What is the corresponding behavior in quantum mechanics? (c) What is 
the quantum analogue of the KAM theorem? (d) The construction of Poincare sections 
and Lyapunov exponents involves long time trajectories. Likewise the quantum analogue 
of Poincare sections and Lyapunov exponents will require long time evolutions. 

In the following we work in imaginary time. In order to avoid notational confusion, we 
continue to use the time parameter T. In particular, T —>■ denotes the high temperature 
limit, where the quantum action approaches its classical counter part. T — > oo denotes 
the zero-temperature limit, where the quantum action is determined by the Feynman-Kac 
formula. In Ref. |27|] the following analytic transformation law at temperature zero has been 



derived in 1-D relating the classical action S = J dt^mx 2 + V(x) to the quantum action 
S = Jdt^mx 2 + V(x), 

2m{V{x) - E gr ) = 2m(V(x) - v ) ~ -^2£2fL^sgn(x) . (3) 

2 y/2m(V(x)-v Q ) 

Moreover, the following relation between the ground state wave function and the quantum 
action holds 

Vv(*) = ^e~ ^ d *' E gr = v , (4) 

Eqs.(|^J^) are valid under the assumption of a parity symmetric potential V(x) with a unique 
minimum at x = 0, V(0) = Vq. Similar laws hold for rotationally symmetric potentials in 
D=2 and 3. For a number of systems, where the ground state energy and wave function is 
known analytically, it has been shown [E7J that the quantum action reproduces exactly the 
ground state energy and wave function. Examples in 1-D are: Monomial potentials and 
the inverse square potential. Examples for radial motion in 3-D are: Monomial potentials 
and Coulomb potential (hydrogen atom) in the sector of non-zero angular momentum. In 
particular, the ground state energy coincides with the minimum of the quantum potential 
V and the location of the minimum of the quantum potential coincides with the location 
of the maximum of the wave function. One should note that the transformation law, 
Eq.(^), determines mass x potential, but not both, mass and potential individually. The 
determination of the renormalized mass requires an non-perturbative numerical calculation. 
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In the following we will try to estimate the precision of numerical calculations of the 
quantum action in the zero-temperature limit (long time limit) by contrasting it with the 
analytic expressions of Eqs.(^j^). We consider in 1-D the classical action 



v 2 x 2 + v±x , 



S = J dt^x 2 + V(x), V(x) 
m — 1 
v 2 = 1 

f 4 = 0.01 . (5) 
The Q.M. Hamiltonian has the following ground state energy, 

E gr = 0.710811 . (6) 
For the quantum action we make an ansatz 

S = J dt—x 2 + V(x), V(x) =v + v 2 x 2 + v 4 x 4 + VqX 6 . (7) 

From a fit of the quantum action to Q.M. transition amplitudes (for more details see 
Ref.||26|) we find the action parameters at T = 4.5 

m = 0.9990(7) 
v = 0.79863(1) 
v 2 = 1.013(2) 
v A = 0.0099(9) 

v Q = 0.0000(3) . (8) 

The parameters rh, v 2 , #4, ^6 have asymptotically stabilized at about T = 4.5. However, 
the parameter v has a slow fall-off behavior for T — > 00. For large T we have fitted Vq to 
the numerical data by 

v (T)^t-,ooA+^ + ^ , (9) 

and find A = 0.710819. Thus the asymptotic extrapolation v (T — > 00) = 0.710819 is very 
close to E gr , given by Eq.(||). This confirms E gr = v , predicted by Eq.(f|), i.e. the ground 
state energy is given by the minimum of the quantum potential. Let us further check the 
consistency of other parameters of the quantum potential. This can be done by using the 
transformation law Eq.([3p. Because this transformation law does not determine rh, we have 
taken rh from the numerial data. The transformation law then determines the parameters 
of the quantum potential. We find, 

v 2 = ^[ m E gr /h} 2 = 1.01151 
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2 Jlmvo r 

v 4 = - z V2 - v 2 —] = 0.009967 

3 n m 

C 6 = -^ + -^-^^4-^4^=2.89 10-'. 10 
4 1> 2 5 n m 

This is consistent with the numerical results, given by Eq. ([§]). Finally, let us look at 
the ground state wave function. As reference solution, we have solved the Schrodinger 
equation with the original Hamiltonian. Secondly, we computed the wave function from 
the quantum action via Eq.(|4]), using the numerical data given by Eq.(S). The comparison 
is given in Fig.[l], showing good agreement. In summary of this section, we have shown 
for the zero-temperature limit (regime of the Feynman-Kac formula) that the transition 
amplitude can be parametrized by the quantum action. A comparison of analytical versus 
numerical computation of the quantum action has shown good agreement. 



4 Quantum chaos in 2-D anharmonic oscillator 

As is well known 1-D conservative systems with a time-independent Hamiltonian are inte- 
grable and do not display classical chaos. An interesting candidate system to display chaos 
is the K-system, corresponding to the potential V = x 2 y 2 . It decribes a 2-D Hamilto- 
nian system, being almost globally chaotic, having small islands of stability [|30|. However, 
according to a suggestion by Pullen and Edmonds |[31|| , it is more convenient from the nu- 
merical point of view to consider the following system, which also exhibits classical chaos. 
It is a 2-D anharmonic oscillator, defined by the following classical action, 

S = f dt \ m (.x 2 + V 2 ) + v ( x iV)i V(x,y) = v 2 (x 2 + y 2 ) + v 2 2X 2 y 2 
Jo I 

m = 1 
v 2 = 0.5 

V 2 2 = 0.05 . (11) 

We work in imaginary time and the convention ft = ks = 1 is used. For the corresponding 
quantum action, we make the following ansatz, compatible with time-reversal symmetry, 
parity conservation and symmetry under exchange x *-> y, 

S = / dt-rh(x 2 + y 2 ) + V(x,y), 
Jo I 

V = v + v 2 (x 2 + y 2 ) + v 2 2X 2 y 2 + v 4 (x 4 + y 4 ) . (12) 

We have determined numerically the parameters of the quantum action for transition times 
from T = up to T = 4, where a regime of asymptotic stability is reached. For T = 4 
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corresponding to temperature r = 0.25, we found 



rh = 0.99814(6) 

v = 1.29373(2) 

v 2 = 0.5154(3) 

v 22 = 0.0497(1) 

v A = -.0008(7) . (13) 

We have included terms xy, xy, xy 3 + x 3 y. Their coefficients were found to be very small 
(compared to machine precision). Also we included terms x 2 y 4 + x 4 y 2 , x 4 y 4 . Those coef- 
ficients were found compatible with zero within error bars. The quantum action slightly 
modifies the parameters v 2 and the parameter v 22 . 

Like it was done in sect . for the 1-D case, we would like to check also in the 2-D 
system how well the quantum action reproduces the ground state properties in the zero- 
temperature limit. The classical action considered now is neither rotationally symmetric, 
nor integrable. Nevertheless, we can express the ground state energy and wave function in 
terms of the quantum action, and compare it with the solution of the Schrodinger equation. 
We use the property of V(x) having a unique minimum at x = 0. In the limit T — > oo using 
conservation of energy, —T kin + V = e = v , the quantum action can be evaluated along 
its classical trajectory, with boundary conditions x d {t = 0) = Xi n and x d {t — T) — Xfi to 
yield 

i-T r T 



t = S[x d ]fjf =J^ dtf kin + V = J o dt2f kin + e 

— >t^oo v T + rh dt (—7— ) 2 = vqT + m [ + 

Jo dt Wo Ji 



T/2 dt J 

v T + m( f° + f Si d~x d ■ &)) . (14) 



lo dt 

Using the conjecture, Eq.([l|), gives the transition amplitude 

G(x fi ,T;x in ,0) — > T ^oo Z e ~ v ° T/h exp[-m f * dx d ■ ^] exp[-m f dx ■ ^] , (15) 

Jo dt Jx in dt 

where we subsumed Z / ' Zq into S and Zq is a constant of dimension 1/L D . Comparison 
with the Feynman-Kac formula yields the ground state energy and wave function, 

E gr = v , 4> gr (x) = ^- exp[-m / dxd ■ -^f] , (16) 



N Jo dt 

where x d is the trajectory from x d {t = T/2) = o to x c [(t — T) — x. The ground state wave 
function thus obtained from the quantum action has been compared to the solution of the 
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Schrodinger equation. This is shown in Fig. [2] . We also compared the ground state energy. 
The solution of the stationary Schrodinger equation yields E gr = 1.01207. Like in the 
1-D system (see Eq.(fJ)) the ground state energy is given by the minimum of the quantum 
potential. Extrapolation of vq to large T yields vq — ► 1.0127. One observes that the 
quantum action describes well the physics at zero temperature also of the 2-D anharmonic 
oscillator. Compared to the 1-D system, however, the error in the wave function is larger 
by one order of magnitude. A reason is certainly due to the fact that in 1-D we have an 
analytical expression @ for the wave function from the quantum action, while in 2-D this 
is not available. 

Now let us look at the chaotic behavior of the classical and the quantum action. We 
have computed Poincare sections from the classical action and from the quantum action for 
a variety of temperatures. The equations of motion of the quantum action have been solved 
using a 4-th order Runge-Kutta algorithm, and Henon's algoritm was used to compute the 
Poincare sections. Let us discuss what we expect. First, the classical system under consid- 
eration is chaotic due to the anharmonic interaction term. This term is small compared to 
the harmonic oscillator term. This property proliferates to the quantum potential. Here 
"smallness" is meant in the following sense: The quantum system at zero temperature is 
determined by the ground state properties. We have computed r = (r) gr , the radius of the 
ground state wave function (corresponding to the Bohr radius ao in the hydrogen atom). 
We obtain tq = 0.8766. It turns out that the anharmonic term is smaller than the harmonic 
term for any length r < r . Thus one expects roughly the same "amount" of chaotic behav- 
ior in the quantum system as in the classical system. Second, renormalisation (i.e. change 
of action parameters when going from the classical action to the quantum action) tends to 
increase harmonic terms, i.e. drive the system to a Gaussian fixed point. At a Gaussian 
fixed point the system is integrable. i.e. non-chaotic. This effect has been observed in a 



number of numerical studies in 1-D systems [26, 28]. In a system with a symmetric double 



well potential this effect was found to lead to a softening of the double well shape of the 
quantum potential (lower barrier, more narrow minima, even vanishing of the double well). 
This translates to a softening or even an "evaporation" of quantum instantons compared 
to the classical instantons . This effect is present also in the system under consideration 
here. We find Sv 2 = v 2 — v 2 = 0.0154(3) compared to 5v 22 = v 22 — v 22 = —.0003(1). We 
expect that this effect may lead to a slight softening of the chaotic behavior in the quan- 
tum system. Third, a weak anharmonic term means that the quantum potential reaches a 
regime of asymptiotic stability as function of T already for small T (T m 4). That is, when 
going from high to low temperature, already at relatively high temperature the asymptotic 
regime (Feynman-Kac limit) is reached. 

We computed the Poincare sections for the quantum action at temperature r = 0.25 
(corresponding to T = 4) for energies E = 10, 20, 50. One should keep in mind that the 
classical action at T = is equivalent to a quantum action at temperature r = oo. Those 
Poincare sections compared to their counter part from the classical action are shown in 
Figs. [3-5]. One observes that the quantum system also displays chaos, and the Poincare 
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sections are slightly different from those of the classical action. Like in the classical case 
also in the quantum system the "amount" of chaoticity increases with increasing energy. 
Moreover, the difference between classical and quantum Poincare sections becomes more 
accentuated with increase of energy. For comparison, we also display in Fig. [6] the Poincare 
sections of the quantum action at E = 50 for temperature r = 2 and r = 0.5, i.e. closer to 
the high temperature limit. 

5 Summary 

We have discussed the use of the quantum action, which can be considered as a renor- 
malized classical action at finite temperature. Based on the quantum action we suggest a 
proper definition of quantum chaos and explore it numerically. As example, we have con- 
sidered in 2-D the harmonic oscillator with a weak anharmonic coupling and computed the 
quantum action at various temperatures. We compared Poincare sections of the quantum 
action with those of the classical action and found some difference in their chaotic behavior. 
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Figure Caption 

Fig.l 1-D anharmonic oscillator. Ground state wave function from Schrodinger equation 
(open dots) and from quantum action (squares). Their difference is shown in the 
lower part. 

Fig. 2 Ground state wave function of 2-D anharmonic oscillator, x indicates the location of 
cuts through the wave function. The upper part shows the solution from Schrodinger 
equation (full line) and from quantum action (different symbols). The difference is 
shown in the lower part. 

Fig. 3 2-D anharmonic oscillator. Poincare sections at energy E = 10. Upper part classical 
action, lower part quantum action at temperature r = 0.25. 

Fig.4 Same as Fig. [3], but E = 20. 

Fig.5 Same as Fig. [3], but E = 50. 

Fig. 6 Poincare sections of quantum action at E = 50, r = 2 and r = 0.5. 
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